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In contrast to the recent prediction of high px hadron suppression 
within the parton saturation model, it is shown that multiple parton 

d. 

scattering suffered by the projectile will enhance high px hadron spectra 
in d + A collisions relative to a superposition of binary p + p collisions 
at RHIC. A stronger enhancement in the forward rapidity region of the 
O ■ projectile is also predicted, resulting in a unique rapidity asymmetry of 

o i 

£^ . the hadron spectra at high px- The shape of the rapidity asymmetry 

, should be reversed for low px hadrons that are dominated by soft and 

coherent interactions which suppress hadron spectra in the projectile 
O , rapidity region. The phenomenon at the LHC energies is shown to be 



qualitatively different because of parton shadowing. 



X 
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I. INTRODUCTION 

Recent experiments [1,2] at the Relativistic Heavy- ion Collider (RHIC) have shown 
a significant suppression of high px hadron spectra in central Au + Au collisions 
that was predicted [3,4] as a consequence of parton energy loss or jet quenching in 
dense matter. In addition, the same mechanism is predicted to produce azimuthal 
anisotropy in high px hadron spectra [5] that was also observed in experiments at 
RHIC [6]. This is a dramatic departure from the heavy- ion collisions at the SPS 
energies where no significant suppression of high px spectra is observed [7,8]. Since 
theoretical studies [9-13] of parton propagation in a dense medium all show that the 
parton energy loss induced by multiple scattering is proportional to the gluon density, 
RHIC data thus indicate an initial gluon density in central Au + Au collisions at the 
RHIC energies that is much higher than that in a large cold nucleus [14]. 

More precise extraction of the parton energy loss from the final hadron suppression 
in A + A collisions, however, requires the understanding of normal nuclear effects in 
p + A collisions. As pointed out by many early [15-17] and recent [18-20] studies, the 
high px hadron spectra can also be modified by initial multiple scatterings in p + A 
and A + A collisions giving rise to the observed Cronin effect. Within a multiple 
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scattering model, the high px hadron spectra are normally enhanced relative to p + p 
collisions, except in the kinematic region where the EMC effect [21] is important (The 
EMC effect is the depletion of parton distributions in x ~ 0.2 — 0.8 in nuclei caused 
by the nuclear binding effect). Such a normal nuclear enhancement will not affect the 
interpretation of the hadron suppression in central A + A collisions as a consequence 
of jet quenching, though inclusion of it is important for more precise extraction of 
parton energy loss and the initial gluon density. 

A very different mechanism was recently proposed for the observed high px hadron 
suppression based on the parton saturation model [22] which also predicts a similar 
suppression in p + A collisions at RHIC, in contrast to the predicted enhancement by 
the multiple parton scattering model. While such a model is not yet checked against 
the existing p + A collisions for energies up to y/s = 40 GeV where enhancement of 
hadron spectra for p T > 2 GeV/c has been successfully explained by the multiple 
scattering model, the up-coming data of d + A collisions at RHIC will attest the 
relevance of parton saturation at the RHIC energies. 

In this letter, we will point out an additional feature in the rapidity dependence 
of the Cronin enhancement due to multiple parton scattering in d + A collisions. 
Such a rapidity dependence was also studied recently by Vitev [23]. However, we 
predict here a unique rapidity asymmetry of high p? hadron spectra due to stronger 
Cronin enhancement in the forward (projectile) region as a result of the transverse 
momentum broadening of the initial partons inside the projectile. The shape of 
the rapidity asymmetry will also depend on the nuclear modification of the parton 
distributions inside a nucleus, in particular at small x. As one decreases px, the parton 
shadowing will reduce the hadron spectra in the forward region, thus changing the 
rapidity asymmetry. When soft and coherent interactions become dominant at very 
low pr < 1 GeV/c, the shape of the rapidity asymmetry will be reversed because 
of the strong suppression of hadron production in the projectile region relative to a 
superposition of binary p + p collisions. We will calculate the rapidity asymmetry 
within a perturbative QCD (pQCD) parton model and study the effect of nuclear 
modification of parton distributions. 
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FIG. 1. Pseudo-rapidity distribution of charged hadrons in minimum-biased d + Au col- 
lisions at y/s = 200 GeV from the HIJING [25] model. 

In a Glauber multiple parton scattering model, large pt spectra are generally en- 
hanced relative to the binary model of hard scattering. It can be shown [24] that 
a combination of absorptive corrections and the power-law behavior of perturbative 
parton cross section leads to a nuclear enhancement at high p T that decreases as l/p T - 
The same absorptive processes suppress the spectra relative to the binary model at 
low p T , where soft processes dominate and the pt spectra deviate from a power- law 
behavior. Since the pQCD model cannot be applied to soft processes, one has to 
resort to phenomenological models like the string model in which coherent particle 
production is modeled by a string excitation for each participant nucleon. Shown in 
Fig. 1 is the rapidity distribution (integrated over transverse momentum) of charged 
hadrons in minimum-biased d + Au collisions at y^s = 200 GeV from the HIJING 
Monte Carlo model [25], which employs the string model for soft and coherent interac- 
tions. Because most of the soft particles are produced through string fragmentation, 
their number should then be proportional to the number of participants. The rapid- 
ity distributions of hadrons from the string should also follow their parent nucleons. 
Since there are more target nucleon participants than projectile nucleons in d + A 
collisions, the hadron rapidity distribution is thus asymmetrical with respect to 77 = 0. 
For high pt hadrons, the underlying processes are hard parton scatterings. The multi- 
plicity should be approximately proportional to the number of binary scatterings and 
the rapidity distributions should be approximately symmetric with respect to i] = 0. 
This is roughly the case for high p T hadrons in HIJING model as shown in Fig. 1. 
One, however, can notice some asymmetric effect in the rapidity distribution of large 
Pt hadrons. This is partially due to the coherence between transverse jets and the 
beam remnants, which is responsible for the asymmetric pedestal effect underlying a 
jet event in p + A collisions. In addition, energy and quark number conservation will 
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also cause some asymmetric effects. These effects are most important in the large 
rapidity region. It is expected that they are small in the central rapidity region. 

Though the HIJING model has incorporated hard processes, it has not included 
transverse momentum broadening of initial partons. For this purpose, we use the 
pQCD model as employed in Ref. [17]. We will use a lowest order (LO) pQCD- 
inspired parton model in which the inclusive particle production cross section in pp 
collisions is given by [26] 

da h r 

- K / dx a dx b d 2 k aT d 2 k bT g p (k aT , Q 2 )g p (k bT , Q 2 ) 



dyd 2 pT 



abed 



fa/ P (x a , Q 2 )f b /p(x b , Q 2 ) h/c ^ ~ ^(ab -> cd), (1) 

nz c at 



where D%, c (z c , Q 2 ) is the fragmentation function of parton c into hadron h as param- 
eterized in Ref. [27] from e + e _ data, and z c is the momentum fraction of a parton 
jet carried by a produced hadron. The K 1.5 (at y/s = 200 GeV) factor is used 
to account for higher order QCD corrections to the jet production cross section. The 
parton distributions f a /N(x, Q 2 ) in a nucleon are given by the MRS D— ' parameteri- 
zation [28]. The initial transverse momentum distribution gN(kx,Q 2 ) is assumed to 
have a Gaussian form, 
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9N{k T ,Q 2 ) = -- I -e- k Tl^ N . ( 2 ) 
■K{k T ) N 

Following [29], we choose a Q-dependent average initial transverse momentum, 

(k 2 T ) N (Q 2 ) = 1.2(GeV 2 ) + 0.2a s (Q 2 )Q 2 , (3) 

which should include both the intrinsic and pQCD radiation-generated transverse mo- 
mentum in this LO calculation. The form of the Q-dependence and the parameters are 
chosen to reproduce the experimental data [17], especially at low energies. Following 
the same approach as in Refs. [29,30], we choose Q 2 to be Q 2 = 2siu/(s 2 + i 2 + u 2 ). 

To take account of multiple initial-state scattering, we assume that the inclusive 
differential cross section for large pt particle production is still given by a single 
hard parton-parton scattering. However, due to multiple parton scattering prior to 
the hard processes, we consider the initial transverse momentum kr of the beam 
partons to be broadened. Assuming that each scattering provide a kx kick which 
also has a Gaussian distribution, we can effectively change the width of the initial hp 
distribution. Then the single inclusive particle cross section in minimum-biased p + A 
collisions is, 

-KJ2 d 2 bt A {b) j dx a dx b d 2 k aT d 2 k bT g A (k a T,Q 2 ,b)g p (k b T,Q 2 ) 



dyd 2 pT 



abed ' 



fa/p(Xa, Q 2 )fb/A(x b , Q 2 , 6) ^ — - ^r(Q& ~> C(Z), (4) 



where is the nuclear thickness function normalized to / d 2 btA{b) = A. We will 
use the Woods-Saxon form of nuclear distribution for tAip) throughout this paper 
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unless specified otherwise. The parton distribution per nucleon inside the nucleus 
(with atomic mass number A and charge number Z) at an impact parameter b, 



f a / A (x,Q 2 ,b) = S a/A (x,b) 



jfa/ P (x,Q 2 ) + (l-j)fa/n(x,Q 2 ) 



(5) 



is assumed to be factorizable into the parton distribution in a nucleon f a /N(x,Q 2 ) 
and the nuclear modification factor S a / A (x,b), for which we take both the new HI- 
JING [31] and EKS [32] parameterizations. The initial parton transverse momentum 
distribution inside a projectile nucleon going through the target nucleus at an impact 
parameter b is still a Gaussian with a broadened width 

(k 2 T ) A (Q 2 ) = (k 2 T ) N (Q 2 ) + S 2 (Q 2 )(u A (b) - 1). (6) 

The broadening is assumed to be proportional to the number of scattering u A (b) the 
projectile suffers inside the nucleus, which is assumed to given by 

3A I 

v A {b) = a NN t A {b) = vnn—^Ji - b 2 /R 2 A (7) 

ZTTJrC^ 

in a hard sphere nuclear distribution, where R A = 1.12A 1 / 3 fm and u^n is the inelas- 
tic nucleon-nucleon cross section. We also assume that kr broadening during each 
nucleon-nucleon collision 5 2 also depends on the hard momentum scale Q = P^ et in 
the parameterized form, 

which is chosen to best fit the existing experimental data in p + A collisions [17] up to 
\fs = 40 GeV. The predictive power of this model lies in the energy and flavor depen- 
dence of the hadron spectra. It is straightforward to also calculate hadron spectra in 
d + A collisions, incorporating the initial hp broadening and parton shadowing. We 
have tried different forms of nuclear distribution for deuteron and find little difference 
in the final results. So we will still use the Woods-Saxon distribution for tdib). In 
heavy nuclear A + A collisions, one can incorporate parton energy loss induced by the 
dense medium through modified fragmentation functions [33]. 

The pQCD model described above has been compared to experimental data for 
pp, pp and pA collisions at various energies [17]. One found that both the intrinsic kp 
and the nuclear broadening are very important to describe the existing data in pp and 
pA collisions, especially at around SPS energies. One can find some more detailed 
description of this pQCD-inspired model in Ref. [17]. We will restrict ourselves to 
the study of the rapidity dependence of the nuclear modification in d + A collisions 
at RHIC and LHC energies in this paper. 
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FIG. 2. Calculated rapidity distributions of charged hadrons with large transverse mo- 
mentum in minimum-biased d+Au collisions at \fs = 200 GeV. The dashed lines correspond 
to a superposition of binary N + N collisions. The dotted lines have nuclear broaden- 
ing of initial parton transverse momentum but without parton shadowing. The solid and 
dot-dashed lines use the EKS [32] and HIJING [31] parameterization of parton shadowing, 
respectively. The spectra have been scaled by the numbers in parentheses. 

Neglecting initial transverse momentum, the initial momentum fractions are re- 
lated to the transverse momentum and rapidities of the final jets by x\ = xt{g Vi + 
e y2 )/2, %2 = xT(e~ yi +e~ y2 )/2, xt = 2Et/ yfs, where Et, yi and y-i are the transverse 
momentum and rapidities of the produced jets, respectively. Large positive rapidities, 
therefore, correspond to large parton fractional momentum x\ from the projectile and 
small momentum fraction Xi from the target. Conversely, negative rapidities corre- 
spond to small X\ and large x%- In our parton model calculation, we assume the 
final hadron rapidity to be the same as that of the fragmenting jets. Shown in Fig. 2 
are the rapidity distributions of the charged hadron spectra for four different values 
of pt- Without nuclear modification of the parton distributions (shown as dotted 
lines), the kx broadening of the projectile partons enhances the particle spectra. The 
enhancement is the strongest in the forward (projectile) region, thus giving rise to a 
rapidity asymmetry. The EKS [32] parameterization has a strong anti-shadowing at 
Xi ~ 0.1 — 0.2 for partons from the nuclear target. Such strong anti-shadowing further 
enhances the rapidity asymmetry (shown as solid lines). However, the HIJING [31] 
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parameterization has mostly shadowing in this region and thus reduces the rapidity 
asymmetry in the spectra (dot-dashed line) which is still visible for px > 3 GeV/c. 
At pt > 10 GeV/c, the target parton distribution at x<i ~ 0.2 — 0.8 is suppressed 
due to the nuclear binding, which is known as the EMC effect [21]. This reduces the 
hadron spectra in the backward (target) rapidity region and thus further enhances 
the rapidity asymmetry. 

As one decreases the transverse momentum so that parton shadowing in the target 
becomes significant, the hadron spectra in the forward region are strongly suppressed. 
It can even overcome the enhancement caused by transverse momentum broadening. 
In this case the rapidity asymmetry is reversed, as shown in Fig. 2 for pt = 1.5 GeV/c. 
This is much like the asymmetric rapidity distribution of soft particles produced via 
soft and coherent interactions as shown in Fig. 1. The validity of the parton model 
at pt = 1-5 GeV/c may be questionable. However, it clearly shows the trend of the 
rapidity asymmetry as one changes the value of pt- 
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Calculated rapidity distributions of charged hadrons with large transverse mo- 



mentum in minimum-biased d + Pb collisions at = 5.5 TeV. 



We also show in Fig. 3 similar rapidity distributions for d + Pb collisions at yfs = 
5.5 TeV where one has access to much larger values of transverse momentum. At such 
high energies and transverse momenta, the parton scattering cross is much flatter in 
Pt than at lower energies. Therefore, the final hadron spectra are less sensitive to the 
transverse momentum broadening due to initial multiple scatterings. The rapidity 
asymmetry caused by the kr broadening is thus very weak, as shown by the dotted 
lines. The nuclear modification of the parton distributions in the target is causing 
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most of the rapidity asymmetry in the spectra at these energies. Here we used only 
EKS [32] parameterization of nuclear modification of the parton distributions. The 
HIJING [31] parameterization does not have any scale dependence, which is very 
important at LHC energies since the shadowing still has significant effects on jet 
production with large transverse momentum. 

To demonstrate the px dependence of the effect of parton shadowing and kx 
broadening, we show in Fig. 4 and 5 the nuclear modification factors defined as the 
ratio of charged hadron spectra in d + A over that in p + p normalized to the averaged 
number of binary nucleon collisions, 

J? ( n d(7 AB/ d y d2 PT , Q v 

R ^y) = 2Ad</dydW w 

At the RHIC energy, the transverse momentum broadening enhances the hadron 
spectra for pt = 2 — 8 GeV/c in both forward and backward rapidity regions. The 
exact enhancement in this pt region depends on the parton shadowing and anti- 
shadowing. The enhancement in the forward rapidity is larger than in the backward 
region as demonstrated by the rapidity asymmetry in Fig. 2. For p? > 8 GeV/c, 
the EMC effect of the nuclear modification of parton distribution starts to suppress 
the hadron spectra. Such a suppression is stronger in the target rapidity region than 
in the projectile region. We also find that the suppression due to the EMC effect is 
stronger for kaons than pions. For pt < 2 GeV/c, parton shadowing suppresses the 
hadron spectra in the case of HIJING parameterization. However, the pQCD model 
might not be valid anymore in this small pr region for a quantitative calculation. The 
nuclear modification factor for d + Pb collisions at the LHC energy shown in Fig. 5 
has much smaller variation and the Cronin enhancement is also smaller. 
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FIG. 4. The nuclear modification factor RdA(PT,y) for d + Au collisions at yfs = 200 
GeV. kx broadening and different parameterization of parton shadowing are included. 

In summary, we have studied the rapidity distribution of hadron spectra in high- 
energy d + A collisions, in particular the rapidity asymmetry caused by multiple 
parton scatterings. The effects of initial multiple parton scatterings are incorporated 
via an impact-parameter dependent nuclear shadowing of parton distributions and 
the broadening of the initial fcy carried by partons before they collide and produce 
high pt hadrons. At low pr parton shadowing suppresses the hadron spectra in the 
projectile rapidity region, giving rise to a rapidity asymmetry much like the soft 
particle production via soft and coherent interactions in a string model. However, as 
one increases pt, hadron production is dominated by hard parton scatterings and the 
rapidity distribution is becoming more symmetric. Within a pQCD model, transverse 
momentum broadening via initial multiple scattering enhances the hadron spectra in 
the projectile region at moderate px = 3 — 8 GeV/c, causing a rapidity asymmetry 
opposite to that of soft hadrons. This is in sharp contrast to the prediction of the 
parton saturation model [22]. The coming data of d + Au collisions from RHIC 
can easily distinguish those two models and will be important to verify whether 
jet quenching due to parton energy loss is truly the underlying mechanism for the 
observed hadron suppression in central A + A collisions. 
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FIG. 5. The nuclear modification factor RdA(PT,y) for d + Pb collisions at -y/i = 5.5 TeV. 

We should also caution that our pQCD model cannot take into account the non- 
perturbative effects that can cause some rapidity asymmetry at very large rapidities. 
Such non-perturbative effect is responsible for the asymmetric pedestal effect under- 
lying a jet event in p + A collisions. Our conclusions are more robust within the 
central rapidity region where such non-perturbative effect and constraints by total 
energy conservation are not yet important. 
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